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[0001] Die vorllegende Erfindung betrifft ein Frequenzanalyseverfahren, bei dem aus einem zeitdiskret, In Form von 
Blocken von Abtastwerten vorliegenden Eingangssignal mittels einer Fensterfunktion uberlappend auteinanderfolgen- 

5 de Blocke von Abtastwerten entnommen werden, woraufhin die so mit den Fensterfunktionen bewerteten BIdcke zur 
ErzeugungeinesSatzes von Koeffizienten Fourier-transformiert werden, nachdem Oberbegriff des Patentanspruchs 1 . 
[0002] Es sind hauptsachlich zwel verschiedene Frequenzanalyseverfahren bekannt, mittels derer ein zeltliches, In 
Form von Abtastwerten vorliegendes Eingangssignal in eine Mehrzahl von Koeffizienten, d. h. in den Spektralbereich 
transformiert werden kann. Diese Verfahren lassen sich in sogenannte Tellbandverfahren (vgl. beispielsweise DE 40 

10 11 577 A1) und Tnansformationsverfahren (vgl. beispielsweise US-51 42656) klassifizleren, wennglelch Transformat- 
ions- und Fllterbankanafysen mathematlsch zunachst eng mitelnander venvandt sInd. 

[0003] Fur die Audiosignalanalyse und/oder -codierung ist es zweckmaBig und vortellhaft, die Frequenzanalyse den 
Eigenschaften des menschlichen Ohres als Schallempfanger anzupassen. Eine wesentliche Erkenntnis der Psycho- 
akustik ist hierbel, da3 das Ohr Schallereignlsse, die in eine besteimmte Bandbreite, die sogenannte Frequenzgrup- 
IS penbrelte fallen, gemelnsam bewertet. Schallerelgnisse auBerhalb eIner solchen Frequenzgruppenbrelte werden ge- 
trennt wahrgenommen. 

[0004] Die Frequenzgruppenbrelte Ist - im Gegensatz zu rein technischen Frequenzanalysen - nicht konstant; mit 
zunehmender Mittenfrequenz wird auch ihre Bandbreite groBer. Der Zusammenhang zwischen Frequenzgruppenbrelte 
und Mittenfrequenz Ist daruberhlnaus, wie In Fig. 2 gezelgt, nlchtlinear 

20 [0005] Dabel venwenden die sogenannten Tellbandverfahren zur Spektralanalyse Filterbanke, die eine spektrale 
Zerlegung entsprechend den Frequenzgruppen vornehmen. Ubiicherweise verwendet man hierzu Kaskadenanord- 
nungen aus Halbbandfiltern. Durch die tortgesetzte Halblerung des Basisbandes mittels einer derartigen Halbbandfil- 
terkaskade kann nur eine grobe Annaherung an die Aufteifung in Frequenz gruppen errelcht werden. 
[0006] Die Zeltauflosung von Bandflltern Ist umgekehrt proportional zu Ihrer jeweiligen Bandbreite; sie wird also mit 

25 zunehmender Bandbreite feiner und Ist damit gut den Eigenschaften des Gehors angepaBt. Halbband-Filterkaskaden 
erzeugen aber, bedingt durch das Prinzip der fortgesetzen Halbterung des Basisbandes, relativ groBe Abweichungen 
der Frequenzaufteilung von der Vorgabe durch die Frequenzgruppenauftellung, die Zwicker fur das Gehor gefunden 
hat. 

[0007] FFT-Analysen konnen zwar bei entsprechend feinem Frequenzraster eine genaue Aufteilung in Frequenz- 
30 gruppen gewahrleisten, die Zeltauflosung istdann jedochzu grob und Insbesondere konstant fur alle Frequenzgruppen, 
was wiederum auch nicht den Eigenschaften des menschlichen Schallrezeptors entsprechen kann, wie die V\fahrnehm- 
barkelt von Vorechostorungen bewelst. 

[0008] Prinzlpiell ist die Zeltauflosung durch die Rekonstruktlon eIner Impulsfolge mit Hilfe eines idealen Tlefpasses 
gegeben. Die aus der Impulsantwort eInes Filters ablesbare Zeltauflosung wird als "resolution span" bezelchnet. Sie 
35 1st in Fig. 1 wiedergegeben. 

[0009] Wie in Fig. 1 gezelgt ist, ist die Erzeugung der auf fg bandbegrenzten Funktlon f (t) durch Superposition von 
mit dem jeweiligen Abtastwert bewerteten sinx/k-Funktlonen a^O) mogiich. Die sinx/x-Funktionen (Spaltfunktionen) 
stellen die Impulsantwort eines Idealen Tlefpasses mit der Grenzfrequenz fg dar. Die Zeitauflosung des idealen Tlef- 
passes ist T 

40 [0010] Das Prinzip wird fur ideale Interpolatoren verwendet. Dazu werden sinx/x-Funktlonen mit dem jeweiligen Abt- 
astwert bewertet und um einen Zeitabschnitt T verschoben, der sich aus dem Abstand von Maximum zum ersten 
Nulldurchgang der sinx/x-Funktion ergibt. Die Superposition allerdieser sinx/x-Funktionen ergibt die Interpolierte Zeit- 
funktion. Der Verlauf der auf diese Weise rekonstruierten kontinuierlichen Zeitfunktion stimmt an den Nulldurchgangen 
der Spaltfunktionen mit den Abtastwerten exakt uberein. Diese Bedlngung ist als Interpolatlonsbedingung bekannt. 

45 [0011] Als Zeltauflosung wird eben dieser Abstand von Maximum zum ersten Nulldurchgang bzw. der Abstand der 
Nulldurchgange der sInx/x-Funktion bezelchnet. Es ist leicht einzusehen, daB die Zeltauflosung gleich der Abtastpe- 
riode ist, wenn als Interpolationsfilter ein Tiefpass mit einer Grenzfrequenz fg gleich der halben Abtastfrequenz f^, also 
ein Allpass gewahit wird. Das Produkt aus Bandbreite und Zeltauflosung betragt In diesem Fall B • T = f^ • 1/2fa = 0,5. 
Bel einer bestimmten Bandbreite kann die damIt gegebene Zeitauflosung nicht unterschritten werden. Dies wird in 

so Analogle zu einem ahnllchen, aus der Physik bekannten Sachverhalt als Heisenberg-Grenze bezelchnet. Teilband- 
und/oder Transfomiationsverfahren zur spektralen Zerlegung finden sich in fast alien heute bekannten Codern. Beide 
sind fur eine gehorgerechte Reprasentation im Frequenz- und Zeitbereich ungeeignet. Speziell muB ein MeBverfahren, 
das beanspucht, die Qualitat heutiger Godlerverfahren beurtellen zu konnen, sicher bereits hinslchtllch der Frequenz- 
analyse genauer arbeiten; erst dann ist die Genaulgkeit, mit der hieruber die MIthorschwellen bestimmt werden hoher, 

ss als in den zu beurteilenden Coder-Decode r-Kombinationen. 

[0012] Die Nachteile herkommlicher, heute fast ausschlieBlich angewendeter Verfahren zur gehorrlchtigen Fre- 
quenzanalyse seien in zusammengefaBter Form verdeutllcht: 
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Teilbandverfahren arbeiten zwar an der Heisenberg-Grenze; die Frequenzgruppenbreiten werden jedoch nur grob 
getroffen. 

Transformationsverfahren sind zwar in der Lage, die Frequenzgruppenbreiten im Prinzip beliebig genau nachzu- 
5 bllden. die Zeitauflosung ist jedoch be! weitem zu grob, daruber hinaus konstant und somit nicht gehoradaquat. 

Hybride Frequenzanalysen nahern sich der geforderten Zeit-/Frequenzaufl6sung des Gehors zwar besser. eine 
exakte Nachbildung ist aber auch hiermit prinzlpiell nicht moglich. 

10 [0013] Der Erfindung liegt die Aufgabe zugrunde, ein Frequenzanalyseverfahren zu schatfen, das die Eigenschaften 
des menschlichen Ohres nicht nur im Frequenz- sondem gleichzeitig auch im Zeitbereich nachbildet. Diese Aufgabe 
wird durch ein Frequenzanalyseverfahren gemaB Anspruch 1 gelost. 

[0014] GemaB der Erfindung werden bei einem Frequenzanalyseverfahren der eingangs genannten Art aus einem 
zeitlich in Form von Abtastwerten vorliegenden Eingangssignal mittels Fensterfunktionen uberlappend aufeinander- 
folgende Blocke von Abtastwerten entnommen, worauf die so mit den Fensterfunktionen bewerteten Blocke zur Er- 
zeugung eines Satzes von Koeffizienten Fourier-transformiert werden, wobei aber jeweils ein solcher Block des Ein- 
gangssignales mit einem Satz von unterschiedlichen. symmetrischen Fensterfunktionen gleicher Lange zur Erzeugung 
eines Satzes von mit den jeweiligen Fensterfunktionen gewichteten Blocken von Abtastwerten bewertet wird. Die Fou- 
rier-Transformierten der Fensterfunktionen besitzen, bedingt durch die unterschiedlichen Fensterfunktionen, auch un- 

20 terschiedliche Bandbreiten. 

[0015] Zu jedem der zeitgleich erzeugten Blocke von gefensterten Abtastwerten Ia3t sich uber eine eigene Fourier- 
Transfomriation mindestens ein Koeffizient berechnen, dem die Bandbreite der Fourier-Transformierten dieser Fen- 
sterfunktion zugeordnet ist, wobei die Koeffizienten so gewahit werden, daB die ihnen zugeordneten Frequenzbander, 
bezogen auf ihre Bandbreite, sich nur geringfugig uberlappen oder beabstandet sind. Mit anderen Worten schlie3en 

2S die den Koeffizienten zugeordneten Frequenzbander im wesentlichen luckenlos und uberdeckungsfrei aneinander an. 
[0016] Einbevorzugtes Ausfuhrungsbeispieldes erfindungsgemaBen Frequenzanalyseverfahrens wirdnachfolgend 
unter Bezugnahme auf die beiliegenden Zeichnungen naher eriautert. Es zeigen: 



Fig. 1 eine Darstellung zur Verdeutlichung. da(3 die Erzeugung einer auf Fg bandbegrenzten Funktion f(t) 

durch Uberlagerung von mit dem jeweiligen Abtastwert bewerteten sinx^-Funktionen moglich ist; 

Fig. 2 den Zusammenhang zwischen der Frequenzgruppe AFg in Abhangigkeit von der Mittenf requenz f^; 

Fig. 3 den Zusammenhang zwischen einer Frequenzachse in Hertz als Ordinate und einer Frequenzachse 

3S in BARK als Abszisse; 

Fig. 4 und 5 Frequenz- und Zeitdarstellungen zur Verdeutlichung des Ahnlichkeitssatzes der Fourier-Transfor- 
mation zur Verdeutlichung, daQ zu einer Verbreiterung der Bandbreite eines Spektralkoeffizienten 
eine Verschmalerung der Fensterfunktion im Zeitbereich notig Ist; 

40 

Fig. 6 ein Beispiel fur eine Fourier-Transfonmlerte einer Fensterfunktion, welche die Energieerhaltung ge- 

wahrleistet; 

Fig. 7 ein Beispiel fur eine Fourier-Transformierte, die eine feinere Zeitauflosung, jedoch nur eine redu- 

45 zierte Frequenzauflosung besitzt; 

Fig. 8 ein FluBdiagramm eines BARK-Generators zur Ermittlung von Fensterfunktionen im Zeitbereich, 

die quenzgruppenbrelte Spektralkoeffizienten erzeugen; 

Fig. 9 eine Darstellung zur Beurteilung eines BARK-Generatorentwurfes; 

Fig. 10 eine Darstellung des iterationsergebnisses des BARK-Generators in Abhangigkeit von der Raster- 

weite; 

ss Fig. 11 eine Darstellung des Quotienten aus der Anzahl unterschiedlicher Frequenzgruppen und der Ge- 

samtzahl der Bander in Abhangigkeit von der Rasterzahl; 

Fig. 12 ein Blockdiagramm einer Vorrichtung zur Durchfuhrung der BARK-Transformation auf der Grund- 
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lage einer 512-Punkte-FFT; 

Fig. 13 ein rekursives Filter erster Ordnung zur Bildung des Goertzel-Algorithmus; 

^ Fig. 14 ein rekursives Filter zweiter Ordnung zur Bildung eines Goertzel-Algorithmus zweiten Grades; 

Fig. 15 eine transponierte Struktur eines Filters zur Bildung eines Goertzel-Algorithmus zweiten Grades; 

Fig. 16 eine bevorzugte Austuhrungsform eines rekursiven Filters zur Bildung des Goertzel-Algorithmus 

10 zweiten Grades als Vorrichtung zur Berechnung eines Koeffizienten der BARK-Transformation; 

Fig. 17 eine Darstellung zum Vergleich des Rechenaufwandes fur einen komplexen Spektralkoeffizienten 

bei Venwendung einer schnellen Fourier-Transformation gegenuberdem Goertzel-Algorithmus. auf- 
getragen uber verschiedene Transformationslangen N; 

IS 

Fig. 18 eine Darstellung der Anzahl GAder Goertzel-Algorithmen, die innerhalb des Rechenzeitbedarfs fur 

eine schnell Fourier-Transformation auf einem Signalprozessor gerechnet werden konnen, aufge- 
tragen uber verschiedene Transformationslangen N; 

^ Fig. 19 eine Darstellung des Frequenzganges als Summe der Amplitudenubertragunqsfunktionen; und 

Fig. 20 eine Darstellung der Betragsubertragungsfunktion der Einzelfilter einer BARK-Transformation. 

[0017] Transformationsverfahren lassen sich offensichtlich mit deutlich weniger Aufwand als Filterbanke reallsieren. 
2S Auf der anderen Seite sind sie aber nur dann mit Erfolg einsetzbar, wenn auf der spektralen Achse (der Frequenzachse 
in Hertz) lediglich aqutdistant angeordnete Reprasentanten benotigt werden. Eine nicht aquidistante Anordnung auf 
der spektralen Achse uber eine Gruppierung herzustellen, fuhrt zu einer Realisierung, die von vome herein nicht ge- 
horrichtig 1st, well sie nicht entlang der Heisenberg-Grenze arbeitet. 

[0018] Der Zusammenhang zwischen einer Frequenzachse in Hertz und einer Frequenzachse in BARK ist aber - 

50 jedenfalls fur niedrige Frequenzen - naherungsweise linear Fig. 2 [Zwicker, E.; Feldtkeller, R.: Das Ohr als Nachrich- 
tenempfanger, 2. Aufl., Hirzel Verlag, Stuttgart 1967] laBt dies bererts erkennen. Deutlicher wird der lineare Zusam- 
menhang, wenn man die Frequenz in Hertz uber der Tonheit in BARK auftragt, wie Zwicker dies in Fig. 3 getan hat 
[Zwicker, E.; Feldtkeller, R.: Das Ohr als Nachrichtenempfanger, 2. Aufl.. Hirzel Verlag, Suttgart 1967]. Dann ist er- 
kennbar, daB der Zusammenhang zwischen einer Frequenzachse f in Hertz (Ordinate) und einer Frequenzachse z in 

35 BARK (Abszisse) zumindest im unteren Frequenzbereich bis etwa 5 BARK linear ist; das bedeutet, das hier die Teil- 
bander gleich breit sind. Im Bereich von 0 Hz bis etwa 500 Hz kann also in guter Naherung von einem linearen Zu- 
sammenhang gesprochen werden. In diesem Bereich laBt sich daherfur jede Mittenfrequenz in etwa die gleiche Fre- 
quenzgruppenbreite von 100 Hz angeben. Hier wird zunachst eine Fourier-Transformation vom Zeitbereich in den 
Frequenzbereich eingesetzt. Sie kann hier ohne die ubiiche Gruppierung von Spektralwerten, d.h. mit einem idealen 

40 Produkt von Bandbreite und Zeitauflosung verwendet werden. Man gewinnt dabei den Vorteil, zumindest die ersten 
f unf Spektralkoeffizienten gemeinsam in einem Ablauf (z.B. in einer schnellen Diskreten Fourier-Transformation) be- 
rechnen zu konnen. Ubiicherweise ergeben sich bei der Abbildung einer Frequenzachse von 0 Hz bis etwa 22 kHz auf 
eine BARK-Achse etwa 25 Spektrallinien. Sie konnen als Reprasentanten f Or die frequenzgruppenbreiten Bander an- 
gesehen werden. Die Berechnung der ersten funf Spektrallinien mittels einer Fourier-Transformation hat mit einer 

45 Schrittweite auf der Frequenzachse zu erfolgen, die der in diesem Bereich annahemd konstanten Frequenzgruppen- 
breite, also etwa 100 Hz, entspricht. Bei einer Abtastfrequenz von beispielsweise 48 kHz mu8 demnach mit einer 
480-Punkte Fourier-Transformation gerechnet werden. Von den 240 Spektrallinien konnen dann schon die ersten 5 
venwendet werden, die restlichen 235 Spektrallinien mussen verworfen werden, da sie zu hoheren Frequenzen hin 
keine Frequenzgruppenbrelte mehr reprasentieren (ihre Bandbreiten sind dafOr zu schmal). Fur den geforderten Ent- 

so wurf entlang der Heisenberg-Grenze sind wegen der Reziprozitat von Bandbreite und Zeitauflosung B • T = 0,5 zur 
Erzielung breiterer Bander entsprechend schmalere Fenster im Zeitbereich fur weitere, getrennt zu berechnende Fou- 
rier-Transformation en zu wahlen; diese ergeben dann - wie gewunscht - auch eine bessere Zeitauflosung. 
[0019] Der Ahnlichkeitssatz der Fourier-Transformation besagt, da3 zu einer "Dehnung" der Bandbreite eines Spek- 
tralkoeffizienten eine "Stauchung" der Fensterf unktion im Zeitbereich notig ist. Zeitauflosung T und Bandbreite B ver- 

ss halten sich in reziproker Weise zueinander. Diese Reziprozitat wird mittels der in Fig. 4 und Fig. 5 veranschaulichten 
Korrespondenzen der Fourier-Transformation gezeigt. Zur Erzielung einer hoheren Bandbreite muB demnach die Fen- 
sterfunktion im Zeitbereich gestaucht werden; seine Zeitauflosung wird damit besser. Der dem erfindungsgemaBen 
Frequenzanalyseverfahren zugrunde liegende Gedanke ist es also, die Ordnung der Fourier-Transformation (die An- 



4 



EP 0 624 866 B1 



zahi der Koeffizienten) beizubehalten. gleichzeitig aberdafur zu sorgen, da3 die hoheren Spektrallinien Reprasentan- 
ten fur breitere Bander darstellen. Auf diese Weise erreicht man auch gleiche Steilheiten fur die Filterflanken. Dies ist 
die Voraussetzung fur ein konstantes Gesamtubertragungsvertialten fur den Betragsf requenzgang, das fur den Einsatz 
in MeBverfahren insbesondere aber auch in Codierverfaliren, gefordert werden muB. 
s [0020] Mit einer BARK-Transformation wird eine gehorgerechte Transformation vom Zeitbereich in den Frequenz- 
bereich bezeichnet, wobei eine aquidistante Darstellung der Spektral koeffizienten auf einer BARK-Achse (und nicht 
auf einer Hertz-Achse) erreicht wird. Die hier entwickelte BARK- Transformation unterscheidet sich von anderen (meist 
nur gehorangepaBten) Abbildungen auf die BARK-Achse dadurch, da3 sie entlang der Heisenberg-Grenze und damit 
gehorrichtig arbeitet. 

10 [0021] Eine Fourier-Transformation, die groOere -6 dB-Bandbreiten lief ert, als es dem Abstand A f der Spektrallinien 
entspricht, ist im herkommlichen Sinne naturlich unbrauchbar (Fig. 6). 

[0022] Normalerweise ist die Bandbreite Bn aller Spektral koeffizienten der Fourier-Transformation durch die Ord- 
nung N und die Abtastfrequenz fg vorgeschrieben: B^ = fg/N. Die Betragsubertragungsfunktionen schneiden sich dann 
in IH(f)l = 0,5, d.h. an den -6 dB-Punkten (Fig. 7). Dennoch konnen von einer nach Fig. 6 modifizierten Fourier-Trans- 

is formation nur die Linien fur die BARK-Transformation venvendet werden, deren Bandbreiten bzw MIttenfrequenzen 
in das von Zwicker fur die Frequenzgruppenbreiten gegebene Schema passen. Durch die - zumindest innerhalb eines 
Rasters - freie Wahl der Bandbreite jedes Koeffizienten erubrigt sich eine Gruppierung von einigen Koeffizienten zu 
einer Frequenzgfuppemitall ihren Nachteilen, wiesie bei herkommlichen gehorangepaBten Transformationsverfahren 
notig ist. Weil das Frequenzrasterfur jede neu zu suchende Linie einer BARK-Transformation beibehalten wird, lassen 

^ sich nun Bander aneinanderreihen, deren Breite in beliebig guter Naherung der angestrebten Abbildung auf die BARK- 
Achse entspricht. Ein Entwurf entlang der Heisenberg-Grenze ist dabei gewahrleistet. Die Staffelung ist deshalb "luk- 
kenlos" moglich, weil Bandbreiten gewahit werden, die genau in das durch die Lange der Fourier-Transformation ge- 
gebene Af-Raster passen. Die Summe aller Amplitudenubertragungsfunklionen bleibt dabei dennoch konstant: die 
Flankensteilheit der bei der Staffelung entstehenden Bandfilter ist nur von der Anzahl der Filterkoeffizienten (der Fen- 

25 sterlange), nicht aber von der Bandbreite abhangig. Diese Anzahl ist aber konstant, sie entspricht normalerweise (d. 
h. ohne Polyphasenaddition) der Lange der Fourier-Transformation. Weil alle Flanken die gleiche Steilheit aufweisen, 
treten unerwunschte Schwankungen im Gesamtubertragungsverhaften (Summenfrequenzgang) hier nicht auf. 
[0023] Die in Form von Abtastwerten vorliegenden Fensterfunktionen mussen nicht notwendigerwerse die gleiche 
Anzahl von Abtastwerten aufweisen, die der Lange der auf die Fensterung folgenden Fourier-Transformation ent- 

30 sprichL Vielmehr kann durch eine Verlangerung der Fensterfunktion uber die Transformationslange hinaus eine Ver- 
besserung der Flankensteilheit der Einzelfilter erreicht werden. Das Additionstheorem der Fourier-Transformation er- 
laubt es, auch langere Sequenzen als durch die Transformationslange gegeben in den Frequenzbereich zu transfor- 
mieren. Hierzu wird eine Anzahl von aufeinanderfolgenden Teilblocken aus dem bereits gefensterten langeren Ein- 
gangsblock von Abtastwerten gebildet, nach dem Prinzip der Superposition addiert und gemeinsam in den Frequenz- 

35 bereich transformiert. Dieses Prinzip ist auch als Polyphasen-Filterung bekannt [Chrochiere, R.E.; Rabiner L: Multirate 
Digital Signal Processing, Prentice-Hall. Englewood Cliffs, New Jersey, 1983]. Durch die Polyphasen-filterung wird 
insbesondere die Zeitauflosung nicht berOhrt, da die dem Filterentwurf zugrunde liegende Spaltfunktion die gleiche 
bleibt. 

[0024] Die Erfindung liefert femer ein als BARK-Generator bezeichnetes Iterationsverfahren, das nunmehr eriautert 
40 wird. Die Staffelung der Bander wird nach Vorgabe einiger Parameter einem Rechenautomaten uberlassen. Diese 
Staffelung liefert fur jede Frequenzgruppe einen Satz von Filterkoeffizienten fur das Zeitf enster einer schnellen Fourier- 
Transformation. Es hat sich als gunstig emlesen, die Staffelung ab der halben Abtastfrequenz (d.h. der hochsten 
vorkommenden Signalfrequenz) zu tieferen Frequenzen hin und nicht umgekehrt vorzunehmen. (Eine Staffelung ab 
der niedrigsten im Signal vorkommenden Frequenz ist naturlich auch moglich). Die Staffelung gewahrleistet nicht von 
45 vorne herein, daB am Ende die untere Grenzfrequenz der untersten Frequenzgruppe mit 0 Hz zusammenfallt. Brauch- 
bar sind aber auch Iteratlonen, die untere Grenzfrequenzen fur das unterste Band bis hinauf zu etwa 50 Hz liefem. 
Ein moglicherweise im Signal enthaltener, normalerweise unhorbarer Gleichanteii wird dadurch in die Messung stark 
gedampft eingehen. Audiosignale reichen nicht bis zu einer Frequenz bis 0 Hertz hinab. Im Verlauf der Iteration werden 
also weitere Bander unter die bereits ermittelten gestaffelt. Die Bandgrenzen der frequenzgruppenbreiten Bander, aber 
50 auch deren MIttenfrequenzen fallen dabei auf ein vorzugebendes Raster. 

[0025] Bei der Suche nach den Koeffizientensatzen fur weitere Fensterfunktionen wird nach dem in Fig. 8 darge- 
stellten FluBdiagramm vorgegangen (CBW = Frequenzgruppenbreite nach Zwicker). 

Die obere Bandgrenze ist bekannt. Sie ergibt sich aus der unteren Bandgrenze des bereits als passend ermittelten 
55 Zeitfensters fur die vorherige, hohere Frequenzgruppe, wobei 6dB-Bandbreiten zugrunde gelegt werden; bei Ite- 

rationsbeginn wird die obere Grenzfrequenz auf die halbe Abtastfrequenz als Startwert gesetzt. 

Die Mittenfrequenz des als nachstes zu suchenden Bandes wird nun solange auf dem Raster verschoben, das 
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10 



sich durch die Abtastfrequenz und die Ordnung der Fourier-Transformation ergibt, bis aus dem Ensemble von 
moglichen Bandbreiten eine gefunden wird, die der Frequenzgruppenbreite von Zwicker am nachsten kommt. 1st 
das verwendete Frequenzraster nur fein genug, so wird auch die Genauigkeit, mit der eine Frequenzgruppenbreite 
getroffen wird, beliebig gut sein. 

Aus der gefundenen Bandbreite ergibt sich der Koeffizientensatz fur eine der Fensterfunktionen und damit die 
Zeitauflosung dieses Fensters gemaB der Heisenberg-Grenze, d.h. die Anforderung im Frequenzbereich bestimmt 
dabei den Verlauf der Fensterfunktion im Zeitbereich. Hier wird das Filterentwurfsverfahren nach Tletze, U.; 
Schenk, Ch.: Halbleiterschaltungstechnik, 9. Auflage, Springer-Verlag, Berlin, 1989, eingesetzt. 

Die untere Grenzfrequenz des gefundenen Bandes bildet die obere Grenzf requenz des als nachstes zu ermitteln- 
den Bandes. Wird im Verlauf der Iteration die Mittenfrequenz negativ, so ist das Abbruchkriterium erreicht. 

[0026] Es ist uberraschend, da3 bei der Untereinanderstaffelung der frequenzgruppenbreiten Bander weit weniger 
IS als 25 unterschiedliche Zeitfenster notig sind. Dies Ist einerseits auf die nahezu konstante Frequenzgruppenbreite zu 
tiefen Frequenzen bis etwa 5 BARK zuruckzufuhren (vgl. Fig. 3); andererseits ist aber auch entscheidend, da3 durch 
die Vorgabe eines Frequenzrasters eben nurdiskrete Mittenfrequenzen bzw. Bandgrenzen erlaubt sind. Wie Fig. 9 
zeigt, konnen bereits bei einer Diskretisierung auf 512 Frequenzen bis zur halben Abtastfrequenz die 25 frequenz- 
gruppenbreiten Bander ausreichend genau mit nur 14 verschiedenen Fensterfunktionen nachgebildet werden: Die 
20 Beurteilung des BARK-Generator-Entwurfs ist sofort graphisch am Bildschirm moglich. Die Abweichung der Bandbrei- 
ten von den Frequenzgruppenbreiten nach Zwicker ist bereits bei einem relativ groben Frequenzraster von 512 dis- 
kreten Frequenzen ausreichend klein (Beispiel in Fig. 9 fur 44,1 kHz Abtastfrequenz). 
[0027] Der BARK-Generator leistet die folgenden Aufgaben: 

25 - Der Entwurf ist prinzipiell fur jede Abtastfrequenz moglich, speziell aber fur die gangigen Abtastfrequenzen 44,1 
kHz, 48 kHz und 32 kHz; eine Abtastfrequenz von 44,1 kHz Ist fur Compact-Disc Muslkmaterlal ublich; Abtastfre- 
quenzen von 32 kHz und 48 kHz findet man bei DAT-Medien (Digital Audio Tape). 

Untereinanderstaffelung von frequenzgruppenbreiten Bandern. Die Genauigkeit, mit der der Entwurf den Fre- 
30 quenzgruppenbreiten von Zwicker entsprlcht, Ia3t sich durch die Vorgabe eines Rasters steuem. 

Zu Versuchszwecken lassen sich auch Bandbreiten fur die Tellfilter einstellen, die Bruchteilen von einem BARK 
entsprechen: dies eroffnet die Moglichkeit, das Tonhohenauflosungsvermogen des menschlichen Ohres von etwa 
1/27 BARK nachzubilden. 

35 

Berechnung und Ausgabe der Satze von Filterkoeffizienten fur die benotigten Fenster unter Berucksichtigung der 
Polyphasenaddition in verschiedenen Ausgabeformaten. Eine der Ausgabemoglichkeiten ist eine Konstantenta- 
belle, die sofort als Quelltext in ein Programm zur BARK-Transformation eingebunden werden kann. Die Bezeich- 
nung BARK-Generator ist damit gerechtfertigt. 

40 

Bandbreite und Zeitauflosung der Fensterfunktionen ergeben (ahnlich wie bei der trivialen Bandfilter-Kaskade) 
immer einen Entwurf entlang der Heisenberg-Grenze. 

Eine visuelle Beurteilung der Iterationsergebnisse ist sofort durch die eingebundene Graphikausgabe moglich. 
45 Eventuelle Modifikationen konnen interaktiv erfolgen. 

[0028] Wie Fig. 10 zeigt, lauft die Iteration nicht fur beliebige Frequenzraster, sondem nur bei einigen ausgezeich- 
neten Rasterweiten erfolgreich ab. Die Iteration kann immer dann als erfolgreich gelten, wenn die untere Grenzfrequenz 
des untersten Bandes im Bereich von etwa 0 bis 50 Hz liegt. Es ist gunstig, einen im Signal moglicherweise enthaltenen 

so Gleichanteil, der ja normalerweise unhorbar ist, auch in die Messung stark gedampft eingehen zu lassen. Der Durch- 
lassbereich der untersten Frequenzgruppe mu3 daher nicht unbedingt bis 0 Hz hinabrelchen. Die Gesamtzahl der 
Bander, in die der Frequenzbereich von z.B. 0 bis 22,05 kHz unterteilt wird, nahert sich mrt zunehmender Rasterzahl 
(d.h. mit feiner werdendem Frequenzraster) der Frequenzgruppenzahl von 25. Diese Anzahl ist fur die in diesem Bei- 
spiel gewahlte Abtastfrequenz richtig. 

55 [0029] Fig. 10 zeigt auch, da(3 die Anzahl unterschiedlicher Fensterfunktionen dabei immer unterhalb der gesamten 
Anzahl von Frequenzgruppen bleibt. Hierin liegt ein wesentlicher Vorteil des BARK-Generators: Es sind - fur alle prak- 
tisch relevanten Falle - weniger Bewertungen des Eingangsvektors des Zeitsignals mit einer Fensterfunktion notig als 
durch die Transformation frequenzgruppenbreite Bander erzeugt werden. Bildet man den Quotienten aus der Anzahl 
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unterschiedlicher Frequenzgruppen und der Anzahl der Bander, so ergibt sich fur eine Rasterzahl von etwa 512 ein 
Minimunn (Fig. 11). 

[0030] Hier arbertet die BARK-Transformation ann okonomischsten: Im Verhaltnis wird die geringste Anzahl unter- 
schiedlicher Bander benotigt. Noch kleinere Rasterzahlen sind zwar moglich, fuhren in diesem Beispiel aber zu Ab- 
s weichungen von der Frequenzgruppenbreite nach Zwicker, die nbht mehr befriedigen. 

Ma3nahmen zur Reduktion des Rechenaufwandes 

[0031] Der Rechenaufwand zur Ermittlung einer vollstandigen BARK-Transfomnation ist noch relativ groB, da fur 

10 jedes unterschiedliche Fenster eine neue Diskrete-Fourier-Transfprmation gerechnet werden muQ. DIese Vorgehens- 
weise geht mit einem meist nur begrenzt zur Verfugung stehenden Potential an Rechenleistung offenslchtlich sehr 
unokonomisch urn, da die nrieisten der berechneten Fourier-Koeffizienten verworfen werden; nur sehr wenige Koeffi- 
zienten (in vielen Fallen sogar nur ein einziger) werden fur die BARK-Transformation weiterverwendet. 
[0032] Wie in Fig. 12 angegeben, kann eine BARK-Transformation, deren Grundlage eine 512-Punkte FFT bildet, 

IS bereits mit nur 11 unterschiedlichen Teilfenstem berechnet werden. Bei der Berechnung der Koeffizienten fur die 25 
frequenzgruppenbreiten Teilfilter uber erne schnelie Fourier-Transformation (FFT) mussen aber 99.56% aller FFT- 
Koeffizienten venworfen werden. nur 25 aus 5632 berechneten Koeffizienten werden tatsachlich verwendet. 
[0033] Fig. 12 verdeutlicht die Verschwendung von Rechenkapazitat in einem Beispiel mit 11 unterschiedlichen Fen- 
stern und einer Rasterung auf 512 diskrete Frequenzen bis zur Abtastfrequenz: Fur die untersten 10 Spektrallinien 

20 kann wegen des In diesem Frequenzbereich naherungsweise linearen Zusammenhangs zwischen Frequenzgruppen- 
breite und Mittenfrequenz noch ein und dasselbe Zeitfenster venwendet werden. Fur die folgenden 3 Bander (11 bis 
13) wird ein neues Zeitfenster und damit eine neue gemeinsame FFT benotigt. Dies gilt auch noch fur die darauffol- 
genden 2 Frequenzgruppen; erst ab der ISten Frequenzgruppe aufwarts muB jeweils ein neues Zeitfenster pro Fre- 
quenzgruppe eingesetzt werden. Bei der ersten FFT werden dadurch 502 von 512 Koeffizienten venworfen, dann 509 

2S von 512, 510 von 512 und schlieBlich 8 mal 511 von 512. Insgesamt werden von 11*512=5632 Koeffizienten nur 25 
verwendet. Es mussen also 5607/25=224 mal mehr komplexe Fourier-Koeffizienten berechnet werden, als dann fur 
die BARK-Transformation gebraucht werden. Wie in Fig. 12 jedenfalls beispielhaft skizziert wird, ergabe die Berech- 
nung der BARK-Transformation uber eine Anzahl von schnellen Fourier-Transformatlonen mit jeweils unterschiedli- 
chen Zeltfenstern einen unverhaltnismaBig hohen Rechenaufwand. wenn nicht eine effizientere Methode zur Berech- 

30 nung einzelner Koeffizienten als die schnelie Fourier-Transformation eingesetzt wurde. 

[0034] Eine triviale Moglichkeit, nur einen Koeffizienten einer Fourier-Transformation zu berechnen, besteht darin, 
die Bestimmungsgleichung fur eine diskrete Fourier-Transformation anzuwenden. Mehr Vorteile bietet aber der Einsatz 
des Goerlzel-Algorithmus zur Berechnung einzelner Linien [Goertzel, G.: An Algorithm for the Evaluation of Finite 
Trigonometric Series, Amer Math. Monthly 65 (1968), S.34-35]. 

35 [0035] Ausgehend von der z-Transformation X(z) einer Sequenz x(n) von Abtastwerten der Lange N 

40 AsO 

ergibt sich die Diskrete Fourier-Transformation (DFT) mit 

4S z = W=(e'^)=e'. 

[0036] Uber k werden diskrete Frequenzen ausgewahlt. 

[0037] Mit der Substitution z" = W'^*^ = e'i^^ kann die DFT auch als Polynom angegeben werden. Beispielsweise 
ergibt sich f Or N = 4 ein Polynom dritter Ordnung 



so 



ss 



X(z) = x(3) . z^ + x(2) . z^ + x(1 ) . z + x(0). 
[0038] Uber das Horner-Schema lost man die Multiplikationen mit Potenzen von z auf 

X(z) = ((x(3) . z + x(2)) . z + x(1 )) . z + x(0). 
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Unter Eintuhrung einer neuen Variablen y(n) kann daraus eine lineare Differenzengleichung angegeben wer- 



5 y(n) = z.y(n-1)+x(N-n). 

[0040] Anfangsbedingung ist hierbei y(0) = 0. Das gesuchte Ergebnis fur die Fourler-Transformierte findet man fur 
n = N 

X(z) = y(N). 

[0041] Die lineare Differenzengleichung kann auch als rekurslves Filter erster Ordnung aufgefaQt werden. Die ein- 
gangs gemachte Substitution z" = W"*« wird hierzu wieder ruckgangig gemacht. 

IS 

y(n) = W*" . y(n-1 ) + x(N-n); y(0) = 0 

[0042] Die damit gegebene Anordnung ist in Fig. 13 wiedergegeben. 

20 [0043] Die Eingangssequenz muB allerdings in umgekehrter Reihenfolge vorliegen. 

[0044] Au3erdem wirkt sich nachteilig auf die Stabilitat und den Reclienaufwand fur das Filter aus, daB die Zustands- 
variable y(n-1 ) uber einen komplexen Koeffizienten W'^ auf den Filtereingang zuruckgef uhrt wird. 
[0045] Die Anzahl der Multiplikationen und Additionen ist zwar die gleiche wie bei der direklen Berechnung der 
diskreten Fourier-Transformation, die Drehfaktoren mussen hier jedoch nicht explizit berechnet oder uber eine Dreh- 

25 faktortabelle ("twiddle table") zur Verfugung gestellt werden; durch die rekursive Bereclinungsweise werden sie viel- 
mehr in jedem Durchlauf, ausgehend von einem Anfangswert, fur einen neu hinzukommenden Abtastwert der Ein- 
gangssequenz x(n) weitergedreht. 

[0046] Die Anzalil der benotigten Multiplikationen laBt sich auf die Halfte reduzieren, wenn man die Euler'sche Glei- 
chung einsetzt 

30 

^ ^-j2Ek/N ^ ^Q5(2iik/N) - j sin(27ck/N) 

und auf eine indizierte Notation ubergeht 

3S 

= [cos(2ick/N) - j sin (2nk/N)] • y„.., + x^.„ . 

40 

wobei 

y„., = [cos(27ik/N) - j sin (27tk/N)] - y„.2 + x^.^^, . 

45 

[0047] Nach einigen weiteren Umformungen erglbt sich schlleBlich fur y^ 

y„ - [2 . cos(27ik/N)] . y„.^ +y„.2 = -x^.^ + [cos(27tk/N) + j sin (2nkm)] • x^+i.^ . 

50 

[0048] Durch Vergleich mit der allgemeinen Differenzengleichung zweiter Ordnung 

y(k+2) + c^ . y(k+1 ) + c^ • y(k) = bg . x(k+2) + b^ . x(k+1 ) + b^ . v(k) 

55 

konnen die Filterkoeffizienten fur ein rekursives Filter zwerten Grades in der ersten kanonischen Grundform angegeben 
werden, das sich wie in Fig. 14 gezeigt, darstellen laBt. 
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= -2 . cos(2ick/N); Cq = 1 
b2 = 0; b, = (2jck/N) + j sin (2iik/N); b^j = 1 

5 

[0049] In dieser Form besitzt die Gleichung reelle Koeffizienten fur die Ruckfuhrung, allerdings einen komplexen 
E ingangskoeffizienten. 

[0050] Das rekurslve Fitter zwelten Grades In der zweiten kanonischen Struklur kann - wie in Fig. 15 gezeigt - noch 
transponiert werden. 

10 [0051] Da die Eingangsfolge rein reell ist, aber in umgekehrterReihenfolgeam Einganganliegt. gilt fur den ungeraden 
Anteil der Eingangsfolge 

Mt) = -h(t) 

IS 

[0052] Der ungerade Anteil hy(t) liefert aber nur den ungeraden Imaginarteil der Fourier-Transformation; speist man 
die Eingangsfolge in nicht umgekehrter Reihenfolge ein, muB der Imaginarteil der Fourier-Transformierten lediglich 
negiert werden. 

[0053] Fig. 16 zeigt die endgOltige Version des Goertzel-Algorithmus. Die Eingangsfolge x(n) liegt jetzt in normaler 
20 Reihenfolge am Eingang an. Der Ausgangswert mu3 erst bestimmt werden, wenn der EIngangsvektor x(n) vollstandig 
abgearbeitet ist. Zur Berechnung einer Spektrallinie ist es vorteilhaft, zunachst alle Werte der Eingangsfolge x(n) in 
das Filter "hineinlaufen" zu lassen und erst dann den Ausgangswert y(n) Ober eine komplexe Multiplikation aus den 
Zustandsvariablen A^ und A2 zu berechnen, 

[0054] Mit diesem wenig bekannten Algorithmus Ist es also moglich, einzelne LInien elner diskreten Fourier-Trans- 
25 formation zu berechnen; die dabei venwendete Struktur Ia3t sich auch als rekursives Filter auffassen. Ein besonderer 
Vorteil gegenuber der FFT ist, da3 hier keine Drehfaktortabelle benotigt wind. Au3erdem treten innerhalb des Filters 
lediglich rein reelle Additionen und Multiplikationen auf. Nur am Filterausgang wird eine komplexe Multiplikation be- 
notigt; zur Berechnung eines BARK-Koeffizienten mul3 diese auQerdem nur einnnal berechnet werden. Der Aufwand 
zur Berechnung eines Koeffizienten Ober den Goertzel-Algorithmus ist in der praktischen Realisierung auf einem mo- 
30 dernen Einchip-Signalprozessor TMS320C30 geringer als der fur eine komplette FFT (Fig. 1 7 zeigt den Rechenauf- 
wand fur verschiedende Transformations Ian gen N). 

[0055] Fig. 18 zeigt dazu, daO innerhalb der Rechenzeit, die eine schnelle Fourier-Transformation benotigt, in vielen 
Fallen sogar einlge Goertzel-Algorlthmen gerechnet werden konnen. Wie aus beiden Fig. 17 und 18 hen^orgeht, ist 
es beispielsweise bei einer Transformationslange von N = 512 gunstiger, bis zu 13 Spektral koeffizienten Ober den 

3S Goertzel-Algorithmus einzein zu berechnen, als Ober eine schnelle Fourier-Transformation. Bei der Anwendung des 
Goertzel-Algorithmus in der BARK-Transformation gewinnt man daher einen ganz betrachtllchen Geschwlndigkeits- 
vorteil. Der Rechenaufwand fur das in Fig. 12 angegebene Beispiel reduziert sich dadurch enorm; in den meisten 
Fallen braucht sogar nur ein Spektralkoeffizient berechnet zu werden. Die Berechnung elner schnellen Fourier-Trans- 
formation Ober eine Koeffizientenanzahf, die nicht einer Potenz von 2 entspricht, ist zwar moglich, fOhrt aber im allge- 

40 meinen zu einem unverhattnisma3lg hoheren Aufwand auf grund der komplizierteren Struktur der dabeieinzusetzenden 
Verfahren. 

[0056] Ein weiterer Vorteil des Goertzel-Algorithmus besteht demgegenuber darin, daB beliebige Frequenzrasterun- 
gen Oberdasselbe Verfahren (d.h. Ober das gleiche Programm) gerechnet werden konnen. Dazu Ist lediglich die Lange 
des Eingangsvektors und der Winkel 27tk/N gemaB der zu berechnenden Linie k neu vorzugeben. Obwohl der Goertzel- 
4S Algorithmus als rekursives Filter am Stabilitatsrand arbeitet, treten dennoch keine Stabilitatsprobleme (z.B. Uberlauf 
des Rechenwerks) auf, da der Algorithmus nach der Berechnung eines jeden Koeffizienten der BARK-Transformation 
wieder neu aufgesetzt wird. 

[0057] FOr die BARK-Transformation ergibt sich somit gegenuber dem in Fig. 12 gegebenen Beispiel ein weitaus 
(um den Faktor 224) geringerer Gesamtaufwand zur Ermittlung des gehorrlchtlgen Kurzzeitspektrums: 25 Goertzel- 
so Algorithmen zur Berechnung der komplexen BARK-Koeffizienten der 25 Teilbander und 11 Polyphasenfenster zur 
Einstellung der 11 verschledenen Bandbreiten. 

[0058] Beide Teilalgorithmen erfordern im Kern ausschlieBlich M u It ip I ikations-Addit ions-Ope rat ion en und sind daher 
f Or eine AusfOhrung auf einem Signalprozessor pradestiniert. Die konsequente Anwendung des Prinzips der Rezipro- 
zltat und die Polyphasen-Fllterung fuhren zu den in Fig. 20 fur eine Abtastfrequenz von f^ = 44,1 kHz angegebenen 
ss Betragsubertragungsfunktionen der Einzelfllter einer BARK-Transformation, aufgetragen Ober einer linearen Fre- 
quenzachse nomilert auf die Abtastfrequenz von f^ = 44,1 kHz. Die BARK-Transformation weist ein konstantes Uber- 
tragungsverhalten im Durchlassberefch. eine ausgezelchnete Sperrdampfung und eine hervorragende Weltabselektl- 
on (keine Nebenmaxima) auf. 
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[0059] In Fig. 19 ist der Frequenzgang der Summe der Amplitudenubertragungsfunktionen angegebender BARK- 
Filterkaskade angegeben. 

Besonderherten der BARK-Transformation 

5 

[0060] Die BARK-Transformation besitzt in der hier vorgestellten Form viele ideale Eigenschaften. Die Nachteile 
herkommlicher Frequenzzerlegungen uber Transformations- und Teilbandverfahren konnen damit uberwunden war- 
den. Die Hauptvorteile der BARK-Transfomnalion sind: 

10 - Gehorrichtige Frequenzauflosung: Die BARK-Transformation folgt dem von Zwicker ermlttelten Zusammenhang 
fur die Frequenzgruppenbreite an der jeweiligen Mittenf requenz. Der Fehler, der durcli die Vorgabe eines Fre- 
quenzrasters bei der iterativen Ermittlung der Grenzfrequenzen der f requenzgruppenbretten Teilbander entsteht, 
kann beliebig klein gehalten werden, wenn nur das Raster fein genug ist. 

IS - Feinstmogliche Zeitauflosung: Die Impulsantwort eines (Band-)Filters kann als Fensterfunktlon fur eine Transfor- 
mation aufgefalBt werden. Die Transformation skoeffizienten stellen im Frequenzbereich Reprasentanten fur die 
(Fourier-)Transtormierten des Fensters im Zeitbereich dan Deshalb mu3 sich be! der Variation der "Bandbreite 
eines Koeffizienten" auch die Auflosung im Zeitbereich andem. Bei der BARK-Transformation ist daher fur jede 
Frequenzgruppe das Produkt aus der Bandbrelte B und der Zeitauflosung T konstant und gleichzeitig minimal: B 

20 . T = 0,5. Diese Grenze, die als Heisenberg-Grenze bezeichnet wird, kann von keiner anderen Frequenzanalyse 

unterschritten werden. Der Filterentwurf entlang der Heisenberg-Hyperbel gewahrleistet, da3 weder die Zeitauf- 
losung noch die Bandbreite der Teilbander durch die BARK-Transformation verschlechtert wird. 

Linearphasige Tellfilter: Aufgrund der Linearphasigkeit der Teilfilter und den gezeigten idealen Eigenschaften der 
25 aus der Iteration resultierenden Betragsubertragungsfunktionen kann die BARK-Transformation in dieser Form 

als bestens geeignet zur Nachbildung der Gehoreigenschaft angesehen werden. Frequenzgruppen zu bilden. Die 
linearphasigen Teilfilter erfordern keine Laufzeitentzerrung durch Allpasse. Da Oberdies die Filterlangen fur alle 
Teilfilter gleich sind, erubrigen sich auch Laufzeitanpassungen durch Verzogerungsglieder. 

30 - Rechteckformige Teilbander: Der iterative Filterentwurf uber die Fenstermethode kann interaktiv an einem Rechner 
erfolgen. Der Entwurf gewahrleistet - trotz der Beschneidung der Spattfunktion auf eine endliche Lange - Betrags- 
ubertragungsfunktionen fur die frequenzgruppenbreiten Teilbander, die sich durch einen konstanten Durchlassbe- 
reich und eine gute Weitabselektion wegen derfehlenden Nebenmaximaauszeichnen. AuQerdem wird ein nahezu 
konstantes Uberlragungsverhaften der gesamten Filterbank erreicht. Zu den ausgezeichneten Eigenschaften, vor 

35 allem der hohen Flankensteilheit bei relativ geringem Realisierungsaufwand gegenuber einzelnen Bandfiltem. 

tragt nicht zuletzt die Polyphasen-Filterung bei. 

Okonomische Realisierung: Bei der Realisierung wird ausgenutzt, da3 die Frequenzgruppenbreiten im unteren 
Frequenzbereich nahezu konstant sind. Es muQ daher mit weit weniger Fensterfunktionen im Zeitbereich multi- 
40 pliziert werden, als Teilbander nach der Transformation entstehen. 

[0061] Die mit der BARK-Transformation gefundene Frequenzanalyse uber nahezu rechteckformige Teilbander bil- 
det den Idealen Ausgangspunkt fur die weitere Entwicklung eines gehorrichtigen Me3verfahrens einerseits und ge- 
horangepaBter Codierverfahren andererseits. Das Transformationsverfahren laQt sich daruber hinaus auch auf andere 
45 Anwendungen ubertragen, bei denen einen Transformation aquidistanter Abtastwerte vom Zeitbereich in den Spek- 
tralberelch benotigt wird, wobei die Spektralkoeffizienten auf der Frequenzachse nicht-aquidistant angeordnet sind 
und dennoch die Heisenberg-Bedingung erfullt blelben soli. 

so Patentanspruche 

1 . Frequenzanalyseverfahren 

bei dem aus einem zeitdiskret in Form von Abtastwerten vorliegenden EIngangssignal mittels Fensterfunktio- 
ss nen uberlappend aufeinanderfolgende Blocke von Abtastwerten entnommen werden, worauf die so mit den 

Fensterfunktionen bewerteten Blocke zur Erzeugung eines Satzes von Koefflzienten Fourier-transformiert 
werden, 
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indem jeweils ein solcher Block des Eingangssignals mit einem Satz von unterschiedlichen, symmetrischen 
Fensterfunktionen gleicher Lange bewertet wird, urn einen Satz von mit den jeweiligen Fensterfunktlonen 
gewichteten Blocken von Abtastwerten zu erzeugen, wobei sich die Fensterfunktionen durch die unterschied- 
lichen Bandbreiten ihrer Fourier-Transformierten voneinander unterscheiden, und 

indem jeder der zeitgteich erzeugten Blocke von get ensterten Abtastwerten eines Satzes einer eigenen Fou- 
rier-Transfomiatbn derart untenworfen wird, 

dadurch gekennzeichnet, daS 

zu jedem der zeitgleich erzeugten Blocke wenigstens ein Koeffizient berechnet wird. dem die Bandbreite der 
Fourier-Transformierten dieser Fensterfunktion zugeordnet ist, und 

daB die Koeffizienten so gewahit werden, daB die ihnen zugeordneten Frequenzbander sich, bezogen auf 
ihre Bandbreite. nur geringfugig uberlappen Oder nur geringfugig voneinander beabstandet sind. 

Verfahren nach Anspruch 1, dadurch gekennzeichnet, 

daB die Fensterfunktionen sinx/x-Funktionen sind, deren Nulldurchgange in Abtastzeitpunkte von Abtastwer- 
ten gelegt sind, so daB sich eine Zeitauflosung entsprechend dem gegenseitigen Abstand der Nulldurchgange 
ergibt, wobei die Zeitauflosung umgekehrt proportional zur doppelten Bandbreite der Fourier-Transformierten 
ist. 

Frequenzanalyseverfahren nach Anspruch 2, dadurch gekennzeichnet, 

daB die Fourier-Transformation eine diskrete Schnelle Fourier-Transformation (FFT) ist, die eine der Anzahl 
von Eingangswerten im Zeitbereich entsprechende Anzahl von Koeffizienten im Frequenzbereich lief ert, deren 
Abstand auf der Frequenzachse der Abtastf requenz geteilt durch die Lange der Transfomnation entspricht und 
deren Bandbreite durch die Fensterfunktion vorgegeben ist, und 

daB aus diesen Koeffizienten nur solche selektiert werden, mit denen zusammen mit zeitparallel aus dem Satz 
aus weiteren Schnellen Fourier-Transformationen auf demselben Block von Abtastwerten mit unterschiedli- 
chen Fensterfunktionen zeitparallel zur Verfugung stehenden Koeffizienten bei anderen Mittenfrequenzen sich 
eine Kaskade von Bandfiftern aufbauen laBt, so daB die sich ergebenden Frequenzbander im wesentlichen 
luckenlos aneinander angrenzen. 

Frequenzanalyseverfahren nach Anspruch 2, dadurch gekennzeichnet. 

daB zur Verbesserung der Selektionseigenschaften der so entstehenden Kaskade von Bandfiltern die Lange 
der unterschiedlichen Fensterfunktionen eines Satzes auch ein Vielfaches der Transformattonslange betragen 
kann, 

wobei vor der Transformation die Vielfachheit von Eingangswerten durch Aufteilung und additive Uberlagerung 
(Superposition) der bereits gefensterten Eingangswerte wieder auf die Transformationslange reduziert wird 
(Polyphasenaddition). 

Frequenzanalyseverfahren nach Anspruch 2, dadurch gekennzeichnet, 

daB die selektierten Koeffizienten aus dem Eingangssignal mittels eines Goertzel-Algorithmus zweiter Ord- 
nung berechnet werden. 

Frequenzanalyseverfahren nach einem der Anspruche 1 bis 5, gekennzeichnet durch 

ein iteratives Verfahren zur Bestimmung der diskreten Mittenfrequenzen und Bandgrenzen zur Festlegung 
der Fensterfunktionen, welches folgende Schritte umfaBt: 

Setzen einer Lauffrequenzvariablen (fmme) e'ne Starlfrequenz (foben). die auf dem durch die Lange 
der Fourier-Transformation gewahlten Frequenzraster liegt; 
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- Verschieben der Lauffrequenzvariablen (fmme) s^f dem Frequenzraster urn einen Schritt (df); 

Ermittlung der sich aufgrund einer gegebenen, nicht linearen Funktton der Mrttenfrequenz ergebenden 
Bandbrerte; und 

Wiederholung der Iteration solange, bis eine Bandgrenze des zu der Lauffrequenzvariablen {^miue) 9eho- 
renden Bandes nnit einer Bandgrenze eines vorher ermittelten Bandes zusammenfallt 

7. Frequenzanalyseverfahren nach Anspruch 6, dadurch gekennzeichnet, 

da3 die gegeberie, nichtlineare Funktion der Zusammenhang zwischen der Mittenfrequenz und der Frequenz- 
gruppenbreite ist. 



Claims 

1. Frequency analysis method 

in which window functions are used to extract overlappingly succeeding blocks of sampled values from an 
input signal present In a time-discrete fashion in the form of sampled values, whereupon the blocks thus 
evaluated using the window functions are subjected to Fourier transformation for the purpose of generating a 
set of coefTicients, 

wherein such a block of the input signal is respectively evaluated using a set of different, symmetrical window 
functions of identical length to generate a set of blocks of sampled values weighted with the respective window 
functions, the window functions differing from one another in such a way that their Fourier transforms have 
different bandwidths, and 

characterized that each of the simultaneously generated blocks of windowed sampled values of a set is subjected 
to a dedicated Fourier transformation in such a way that 

for each of the simultaneously generated blocks at least one coefficient is calculated which is assigned the 
bandwidth of the Fourier transforms of this window function, and 

that the coefficients are chosen such that with reference to their bandwidth the frequency bands bands as- 
signed to them overlap one another only slightly or are separated from one another only slightly 

2. Method according to Claim 1 , characterized in that 

the window functions are sinx/x functions whose zero crossings are situated at sampling instants of sampled 
values, resulting in a temporal resolution corresponding to the mutual spacing of the zero crossings, the tem- 
poral resolution being inversely proportional to twice the bandwidth of the Fourier transfomns. 

3. Frequency analysis method according to Claim 2, characterized in that 

the Fourier transformation is a discrete fast Fourier transformation (FFT) which delivers in the frequency range 
a number of coefficients corresponding to the number of input values in the temporal range and whose distance 
on the frequency axis corresponds to the sampling frequency divided by the length of the transformation and 
whose bandwidth is prescribed by the window function, and 

in that of these only those coefficients are selected by means of which it is possible, in conjunction with coef- 
ficients available in a time-parallel fashion from the set of further fast Fourier transformations on the same 
block of sampled values using different window functions, to construct a cascade of band filters in the case of 
other centre frequencies, so that the resulting frequency bands adjoin one another essentially without gaps. 

4. Frequency analysis method according to Claim 2, characterized in that 

in order to improve the selection properties of the cascade of band filters thus produced, the length of the 
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different window functions of a set can also be a multiple of the transformation length, 

It being the case that before the transformation the multiplicity of input values is reduced again to the trans- 
formation length by division and additive superposition of the already windowed input values (polyphase ad- 
dition). 

5. Frequency analysis method according to Claim 2, characterized in that 
the selected coefficients are calculated from the input signal by means of a second order Goertzel algorithm. 

6. Frequency analysis method according to one of Claims 1 to 5, characterized by 

an iterative method for determining the discrete centre frequencies and band limits for the purpose of defining 
the window functions, which method comprises the following steps: 

- setting a running frequency variable (fcentre) a start frequency (f^p) which lies on the frequency raster 
chosen by the length of the Fourier transformation; 

displacing the running frequency variable (fcentre) ^^^P (^) on the frequency raster; 

determining the bandwidth resulting from a given, nonlinear function of the centre frequency; and 

repeating the iteration until a band limit of the band belonging to the running frequency variable (fcentre) 
coincides with a band limit of a previously determined band. 

7. Frequency analysis method according to Claim 6, characterized in that 
the given, nonlinear function is the relationship between the centre frequency and the frequency group width. 

Revendications 

1. Procede d'analyse de frequences 

3S t dans lequel des blocs, qui se suivent en se recouvrant, de valeurs detect^es sont pr^lev^s au moyen de 

fonctions de fendtre k partir d'un signal d'entrSe present sous forme de valeurs d^tectdes de fagon discrete 
dans le temps, k la suite de quoi les blocs ainsi evalues au moyen des fonctions de fenStre sont soumis h une 
transformation de Fourier pour engendrer un jeu de coefficients. 
♦ dans lequel un tel bloc respectif du signal d'entree est lvalue au moyen d'un ensemble de fonctions symetri- 
40 ques differentes de fenetre de meme longueur afin d'engendrer un jeu de blocs de valeurs detect§es ponddres 

au moyen des fonctions respectives de fen§tre, les fonctions de fen§tre se diff^renclant Tune de I'autre par 
les largeurs differentes de bandes de leurs transformees de Fourier, et 
4 dans lequel chacun des blocs engendrSs simultan^ment de valeurs detectees fenetrees d'un ensemble est 
soumis a une transformation de Fourier propre, 



2S 



30 



45 



ss 



caract^ris^ en ce que 



♦ au moins un coefficient auquel est associ§ la largeur de bande de la transform6e de Fourier de cette fonction 
de fenetre est calculi pour chacun des blocs engendr^s simultan^ment, et en ce que 
^0 ♦ les coefficients sont selectionn6s d'une mani^re telle que les bandes de frequence qui leurs sont associ6es 

ne se recouvrent que leg^rement ou sont leg^rement espac6es les unes des autres, en ce qui concerne leur 
largeur de bande. 



2. Proc6d6 selon la revendicatbn 1 , caractSris^ en ce que 

♦ les fonctions de fenetre sont des fonctions sinx/x dont les passages par zero sont places k des instants de 
detection de valeurs ddtect^es, de sorte qu'il en resulte une resolution temporelle qui correspond k la distance 
reciproque des passages par zSro, la resolution temporelle etant inversement proportionnelle au double de 
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la largeur de bande de ia transformee de Fourier. 
Proc§6§ d'analyse de frequence selon la revendication 2, caract6ris6 

♦ en ce que la transformation de Fourier est une transformation rapide de Fourier (FFT) discrete qui fournit dans 
le domaine des frequences un nombre de coefficients qui correspond au nombre des valeurs d'entr^e dans 
le domaine des temps, coefficients dont la distance sur I'axe des frequences correspond k la frequence de 
detection divisee par la longueur de la transformation et dont la largeur de bande est pr^definie par la f onction 
de fenetre, et 

♦ en ce que parmi ces coefficients, seuls sont seiectionn^s ceux au moyen desquels peut §tre constitute k 
d*autres frequences moyennes une cascade de filtres de bande k coefficients disponibles parallelement dans 
le temps au moyen d'autres transformations rapides de Fourier du meme ensemble sur le meme bloc de 
valeurs detectees au moyen de fonctions differentes de fenetres parallelement dans le temps, de sorte que 
les bandes de frequence qui en resultent sont adjacentes I'une k Tautre sensiblement sans lacunes. 

Procede d'analyse de frequences selon la revendication 2, caracterise en ce que 

♦ la longueur des differentes fonctions de fenetre d'un ensemble peut egalement etre egale a un multiple de la 
longueur de transformation afin d'ameiiorer les proprietes de selection de la cascade ainsi produrte de filtres 
de bande, 

♦ la serie des valeurs d'entree etant de nouveau reduite avant transformation a la longueur de transformation 
(addition polyphasee) par division et superposition additive (superposition) des valeurs d'entree dej^ fene- 
trees. 

Procede d'analyse de frequence selon la revendication 2, caracterise en ce que 

i les coefficients seiectionnes dans le signal d'entree sont calcuies au moyen d'un algorithme de Goertzel de 
deuxieme ordre. 

Procede d'analyse de frequence selon I'une des revendications 1^5, caracterise par 

♦ un procede iteratif destine k determiner les frequences moyennes discretes et limltes de bandes pour fixer 
les fonctions de fenetre et comprenant les etapes consistant k : 

• etablir une variable de frequence courante (f^me) ^ frequence de depart (foben) ®st situee sur le 
reseau de frequences seiectionne par la longueur de la transformation de Fourier; 

• decaler d'une etape (df) la frequence variable courante (fmitte) '® reseau de frequences ; 

• determiner la largeur de bande resultante sur la base d'une fonction donnee, non lineaire, de la frequence 
moyenne; et 

• repeter iteration jusqu'^ ce qu'une limite de bande de la bande associee k la variable de frequence 
courante (f^itte) coincide avec une limite de bande d'une bande determlnee precedemment. 

Procede d'analyse de frequences selon la revendication 6, caracterise en ce que 

f la fonction donnees, non lineaire, est la relation entre la frequence moyenne et la largeur du groupe de fre- 
quences. 
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Vorgaben: 
Abtastfrequenz: fa 
Frequenzraster; #lines 




foben = fa/2 
df sfoben/#lmes 




funten = fmitte - (foben - fmitte) 



I 



BW = (foben -funten) | 




Ausgabe an Filterentwurf: 
BW(fmitte) 




foben = funten 
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